arXiv:1509.04164vl [hep-ph] 14 Sep 2015 


A simple merging algorithm for jet angular correlation studies 


Junya Nakamur£[^ 

Institute for Theoretical Physics, University of Tubingen, 
Auf der Morgenstelle I 4 , 72076 Tubingen, Germany. 


Abstract 

A tree level merging algorithm which guarantees the leading order (LO) accuracy of angular correla¬ 
tions between jets is proposed and studied. The algorithm is designed so that n-jet events are generated 
exclusively according to the LO n-parton production cross section and each of the n-jet is close to each 
of the n-parton in terms of the jet measure. As a result, the LO accuracy of angular correlations between 
the n-jet is robust. Furthermore, as long as the n-jet events are exclusively subjects to a study, only the 
LO n-parton production cross section is needed and hence event generation is efficient. Correlations in 
the azimuthal angle difference between the two highest transverse momentum jets with large rapidity 
separations in the top quark pair production are evaluated as examples. The algorithm is validated by 
discussing numerical differences between its predictions and the predictions of a well-established merging 
algorithm. 


^ E-mail: junya.nakamura@uni-tuebingen.de 


1 



Contents 


1 Introduction 

2 Formalism [3] 

3 Algorithm 0 

4 Numerical studies 1101 

5 Summary and discussion 1131 


1 Introduction 

Angular correlations between jets produced together with heavy particles have been studied actively for a 
long time, because they can provide important information about the heavy particles HI HI I IB H HI- 
For instance it has been shown that the distribution of the azimuthal angle difference A(j) = (f>i — (j )2 
between two partons in the gluon fusion production of a Higgs boson plus the two partons is very sensitive 
to a charge-conjugation and parity (CP) property of the Higgs boson [B H H 111 IS]- By observing the 
A(j) distribution and comparing it with theoretical predictions, we can measure CP violation in the Higgs 

sector mm- 

In order to read the information of heavy particles from angular correlations between jets produced in as¬ 
sociation with them, it will be necessary to produce the accurate predictions of observables, such as Acj), which 
measure the angular correlations. Tree level merging algorithms |511^ ITUl[TTl IT^IT^ ITUIT?1 [TO [T71IT51ITO1 [^ . 
which combine leading order (LO) cross sections for multi-parton in final state with the parton shower, 
are nowadays standard tools used for simulating processes including multi-jet in final state. Models of the 
parton shower base the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation 122 US] 
and thus the parton shower guarantees the leading logarithmic (LL) accuracy for the kinematics of produced 
partons. Therefore the accuracy guaranteed in merging algorithms can be seen as the LO plus LL. The 
virtue of merging algorithms is that they can combine LO cross sections smoothly with the LL parton 
shower so that the dependence on an artificial scale at which they are combined is minimized. 

A LO multi-parton production cross section predicts angular correlations between the produced partons 
at the LO accuracy, while the LL parton shower does not have ability to predict angular correlations 
between any partons. Considering this fact, when our objective is to predict angular correlations between 
constructed jets, the accuracy minimally required for the kinematics of the jets should be the LO, neither 
LO-I-LL nor LL. If the kinematics of a jet is determined or largely influenced by the LL parton shower 
during a merging procedure, it hardly has the LO accuracy and hence it is not appropriate to use the event 
containing this jet. Merging algorithms in the literature potentially have the ambiguity in the accuracy of 
jets, namely it is not necessarily clear whether the kinematics of a jet constructed by clustering particles in 
final state after a merging procedure has the LO accuracy or not. This is because their virtue is smooth 
combination of LO cross sections and the LL parton shower. 

In ref |24j . the azimuthal angle difference A(j) = (f>i — (j )2 between the two highest transverse momentum 
prp jets with a large rapidity separation in the ti pair production is studied by using the CKKW-L 
merging algorithm |B [12 [20] with the parton shower model |2B [2B [12 in PYTHIAS [3B [20] ■ There, 
it is found that the correlation between the 2-jet can be lost in a non-negligible fraction of the events 
when the LO cross sections for the ti plus up to 2-parton are merged with the parton shower, because 
a jet originating from the parton shower has a higher pj, than one of the two jets originating from the 
2-parton of the LO cross section. This shows one example of the case that the kinematics of a jet is 
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determined or largely influenced by the LL parton shower during a merging procedure. Although it is 
also found that the loss of the correlation can be avoided by merging the LO ti + 3-parton cross section 
additionally, calculating LO cross sections for higher multiplicity is time-consuming, thus we want to avoid it. 

In this work I construct a new merging algorithm which does guarantee the LO accuracy of angular 
correlations between jets and hence does not have the above ambiguity. The Sudakov suppression is calculated 
in the same way as the MLM [13 H] and the fc_i_-jet MLM |14L I16j algorithms. The differences from these 
algorithms are 

• The definition of jets used during a merging procedure is set identical to the one used during physics 
analyses of jets. 

• The LO n^na^x'P^'^ton production cross section, where denotes the maximal number of partons 

produced by the LO cross section, is not allowed to produce the events which contain more than 

The first point indicates that n-jet events are generated exclusively according to the LO n-parton produc¬ 
tion cross section and furthermore each of the n-jet is close to each of the n-parton in terms of the jet 
measure. As a result, angular correlations between the n-jet should follow those between the n-parton and 
hence the LO accuracy of them should be robust. In addition to this, when n-jet events are exclusively 
subjects to a study, only the LO n-parton production cross section is needed, thus event generation is efficient. 

The second point makes it forbidden to produce the events which contain more than If our 

objective is to predict angular correlations between two jets such as A;/) at the LO accuracy, not only do 
the two jets have angular correlations between them at the LO accuracy, but all of the additional jets also 
should have angular correlations with the two jets at the same accuracy, because the additional jets can 
affect kinematically. If the kinematics of the additional jets is determined by the LL parton shower, 
those jets do not have angular correlations with any other jets. In such a case, the accuracy of A^ will be 
less than the LO. Our objective is to predict angular correlations between jets at the LO accuracy, therefore 
the events which contain more than should not be used for analyses. In my algorithm, this is 

naturally achieved by the above second point. 

The tt production in proton-proton collisions is simulated and the azimuthal angle difference A(j) = 4‘i~4’2 
between the two highest pj, jets with large rapidity separations is studied by using the 2-jet events exclu¬ 
sively. The new algorithm is validated by comparing the A^ distributions with the predictions of one of 
well-established algorithms, the CKKW-L algorithm lidlllD]. It is observed that the new algorithm at 
^max = 2 produces a consistent result with the CKKW-L algorithm at u-max = 3. 

In Section my formalism of tree level merging algorithms and my notations are introduced. Following 
these, in Sectionj^the new merging algorithm and an event generation procedure according to it are described 
in detail. In Section]^ the results of the simulation are presented. In SectionI summarize my hndings. 


2 Formalism 

In this section, my formalism of tree level merging algorithms and my notations are introduced. Following 
these, the new algorithm is described in Section]^ This section is largely based on ref. [24] . 

The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation [^ [23113 be written 
with the Sudakov form factor in the following form |5D] 

^ d q{x,t) _ i , q{x/z,t) 

dt A{t) Jo z ’ A(t) ’ ^ ^ 
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where A{t) denotes the Sudakov form factor 


A(t) = exp ( — 


r* dt' 

,,2 t' In 




( 2 . 2 ) 


Here only the quark parton distribution function q{x,t) (PDF) and the splitting function Pqq{z) for q ^ qg 
without the virtual correction are introdnced in order to simplify writing. The generalization is, however, 
simple. By integrating eq. (2.1) over t it follows that 


q[x,tx) ^ qix,tj^) 

^i^x) ^(^a) 

where a short hand notation is introduced 


dt 




a^dz ~ 

dPqqiz) — 


(2.3) 


(2.4) 


Using eq. (|2.3|) iteratively and dividing it by q{x, t^)/A{tx), we can obtain the following form of the DGLAP 
equation 


qix,tj,) Aitx) 

-hi 


Mx/zi,tj^) A{tx) 
q{x,tx) A{t^) 


+ 


dt 


r-elb) 


dtn 


re(*2) 


/ dp„Jz^) / — 


JO 


dPnq{z 2 ) 


2 JO 


q{x/{z^Z 2 ),tt^) A{tx) 
q{x,tx) A{tj^) 


+ 


(2.5) 


The first term of the right hand side (RHS) in the above equation can be seen as the probability of 
generating no radiation from a quark q{x) in a proton during the scale evolution of the proton between 
tx and {tx > t^). The second term of the RHS represents the integrated probability of generating 
exclusively one radiation from the quark during the evolution, and so on. The left hand side (LHS) ensures 
the probability conservation. 

I write the no radiation probability from a quark q{x) in a proton during the scale evolution of the proton 
between and t 2 (ti > ^ 2 ) the following form [SD] 


^q{ti,t 2 ix) = 


q{x,t 2 ) A(ti) 

q{x,tO A{t 2 )' 


( 2 . 6 ) 


Then, eq. (2.5) can be expressed as 

pt y 


1 —H (t^, a;) + / 


1 JO 


dPqqizO Tlgitx,t^\x) ^qitl,tA\x/z^) 


dti 


e(*i) 


t 


dPqqiZi) 


1 JO 


dt 2 

t. 


2 JO 


q{x,t^) 

#q9(^2)nq(tx: ti; ta; x/z^) 


q{x,t-^) 


q{xl{z^Z2),t2) , , 

qix/z„t2) ^^q[*2.tA,x/(Z,Z2)) 


(2.7) 


It is an easy task to find the explicit form of IlJt^,t 2 ', x) from the above equation |31j 


/ /"‘i dt P'd'i 

n,j(ti,t 2 ;a;) = exp —J dp^^iz) 


,qix/z,t) 

qix,t) 


( 2 . 8 ) 


Given the scale tx and the energy fraction x of an incoming quark in a proton, eq. (2.7) allows us to 
generate radiations from the incoming quark by evolving the proton from the scale tx- This is known as 
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backward evolution [sllEaiao]. 


Eq. (2.7) derived from the DGLAP equation concerns only radiation from an incoming quark in a 
proton i.e. initial state radiation. Here I generalize the equation to the one which can predict radiation 
from outgoing partons i.e. hnal state radiation, as well as initial state radiation. What we should notice 
for this purpose is that the PDFs play a role in constraining scale evolution of the proton, or in other 
words constraining radiation from the quark in the proton during scale evolution of the proton. Hence, 
when radiation from outgoing partons is concerned, we replace the PDFs with a function which bases 
the kinematic information of the outgoing partons. The function constrains radiation from the outgoing 
partons, through the energy and momentum conservation for instance. 


I let {p}x+n denotes a complete specification of an event sample consisting of X + n partons The 
information of two incoming partons is implicitly included. Then, I introduce a function for the evolution of 
a {P}x+n 


f{z,t-Ap}x+n), 


(2.9) 


which constrains the evolution of the {p}x+n at the evolution scale t and the energy fraction z. By using 
the constraint function, eq. (2.7) can be generalized to 

i=ii{tx,t^;{p}x) + / -r ^pAi) /(^i>ti;Mx) ^{h,tAAp}x+i) 

£i Jo 

dt^ A 


t 


dp{zAfd{tx,h] Mx) /(^i: ti; {p}x) 


and accordingly 


1 JO 
/■‘i (jii 

' / dp(z2)n(ti,t2;Mx+i)/(^2>t2;Mx+i) n(t2,tA;Mx+2) 

^2 Jo 


n(ti,t 2 ;Mx+n) = exp dp{z)f{z,t;{p}x+n) 1 ^ 


( 2 . 10 ) 


( 2 . 11 ) 


which is defined as the no radiation probability for a {p}x+n as a whole, during the scale evolution of it 
between and t 2 (ti > ^ 2 )- ta other words, this is the probability that a {p}x+n remains the same during 
the evolution. The splitting probability has the following form 


dp{z) = ^^P{z) 

Z ZTT 

dpiz) = dz^p{z) 


for initial state radiation, 
for final state radiation. 


( 2 . 12 ) 


where the appropriate splitting function(s) should be used for P{z) according to a branching process 
{p}x+n {p}x+{n+i)- initial state radiation, the constraint function f{z,t;{p}x^n) always includes 
the PDFs. The soft gluon singularity at 0 = 1 in the splitting functions will be avoided by introducing 9 
functions in the constraint function. 


Let us consider a hard process which produces X and express its leading order (LO) cross section as 
^This expression is inspired by refs. dsiiis]. 
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cr{X). The DGLAP evolution of the hard process can be expressed by multiplying a{X) by eq. (2.10) 


cr(A) = a{X) U{tx,tj^-{p}x) 

r^x fii /-i 

+ <^ix) dp{z^)n{tx,h\{p}x) ^(ti^txAp}x+i) 

JQ 

f^x /*! 

+ / 7 ^ / dp{zi)^{tx,h'Ap}x) fizi,h'Ap}x) 

Jo 

r^\ cii /*^ 

x/ dp{z2)Il{tj^,t2\{p}x+l)f{z2,t2',{p}x + l)'^{t2,tA:{p}x+2) 

Jtj^ ^2 Jo 

+ •••• (2.13) 

Tree level merging algorithms are designed as to improve the above DGLAP evolution by replacing the terms 
constructed by the leading order cross section times the universal radiation probability with the exact LO 
cross sections for multi-parton in final state [8]. In my notations, it proceeds as 

cr(X) = cr(A) n{tx,tj^-{p}x) 

+ a[x+ i-,{p}x+i >tj^ ,G ;f {zi,ti^-,{p}x) n(G,{rfjc+i) 

+ (j{x + 2 \{p}x+ 2> f'{^i^^i\{p}x) n(ti,t2;{rix+i) 

X f'{z2,t2]{p}x + l) n(t2i {p}x+2) 

+ •••, (2.14) 


The soft and collinear divergences in LO cross sections are regularized in the following way. By using the 
definition of the evolution variable t, calculate the minimum value from a {p}x+n require > G ■ 
This is expressed as {p}x+n > the above equation H Notice that the constraint functions in eq. (2.14) 
are different from those in eq. (2.13), since some part ofthe constraint is already included in the exact LO 
cross section. I call /'(z„_|_]^, {p}x+n) ^ weight function hereafter, since it is used to re-weight an event 

sample during a merging procedure. The cut off scale is called merging scale and separates phase space 
into two regions. Partons in the phase space above the merging scale are produced following exact LO cross 
sections and those in the phase space below the merging scale are produced following the DGLAP equation. 
The further evolution of each term in the above equation will be performed between the scale and a 
lower cutoff scale for complete event generation. In my notations, the term consisting of the cross section 
a{X + n) in the RHS of the above equation is called the n th term. 


Below tree level merging algorithms which have been proposed and studied in the literature are briefly 
reviewed. This information can be useful when I describe my algorithm in Sectionj^ There are fundamentally 
two different types of merging algorithms ROne of the two is constructed to generate events according to the 
improved DGLAP equation strictly, eq. (2.14) in my notations, by explicitly calculating the Sudakov form 
factors. This procedure is necessary in order to minimize the dependence on the merging scale and smoothly 
combine LO cross sections with the parton shower. Algorithms of this type include the GKKW with the 
vetoed parton shower um, the GKKW with the truncated shower [HI[TH] and the GKKW-L [9l[ElEO] 
algorithms. The other of the two types does not calculate the Sudakov form factors explicitly and hence 
event generation follow the improved DGLAP equation approximately. Algorithms of this type include the 
MLM |13[I14] . the fcj_-jet MLM |141I16) and the shower kj_ US] algorithms. In the MLM and the A:j_-jet MLM 
algorithms, partons in final state obtained after the complete shower evolution are clustered according to the 
cone jet algorithm and the exclusive A:_L-jet algorithm, respectively. Then, if the number of the constructed 
jets is not equal to the number of partons produced by a LO cross section or if the constructed jets are not 
close to the partons in terms of the jet measure, the event is vetoed. In the shower kj_ algorithm, if the scale 
of the first emission is above the merging scale, then the event is vetoed. 


^Here it is assumed that the hard process cross section o-(X) is finite everywhere in its phase space. 

^In this work I concentrate only on discussing a type of merging algorithms which introduces the merging scale. However, 
it would be also interesting to study jet angular correlations with other types of algorithms such as m- 
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3 Algorithm 

In this section, the ideas of the new merging algorithm and an event generation procedure according to the 
algorithm are described in detail. 

The new algorithm is constructed with a goal of eliminating the ambiguity in the accuracy of jets which 
potentially exists in merging algorithms as discussed in Section and hence of guaranteeing the leading 
order (LO) accuracy of angular correlations between jets. The following idea is strictly implemented in the 
algorithm: 

The n-jet events are generated exclusively according to the LO n-parton produetion cross section and 
furthermore each of the n-jet is close to each of the n-parton in terms of the jet measure. 

For an algorithm which implements the above idea, which I call the model A algorithm eq. 
be written in the following form: 

0'(-^)ModelA 

= (t{x) n(tj^, {p}x —>■ {rix+ojet; {p}x) 

+ a(^X + 1; {p}x+i > {p}x+l {p}x + liet'Ap}x+l) 

+ a^X + 2 ; {p}x+2 > Qraf) n(tjc_l_2, Mx+ 2 ^ MA+2jets; Mx+ 2 ) 

H- 

+ a(^X + n; {p}x+n > Q™?) ^{'^X+n^ {p}x+n {rfx+njets; {p}X+n) 

H- 

+ cr^x + > Qcut) ^{^X+n^^^Ap}x+n^^^ ^ • (3-1) 

Here the weight functions are omitted in order to simplify writing. In the numerical studies in Section 
they are taken into account. The expression {p}x+n > indicates that the soft and collinear divergences 
in LO cross sections are regularized by the definition of The definition of must respect the 

definition of jets used in analyses, namely a jet clustering algorithm and parameters that the algorithm 
contains. Throughout of this work, the anti-Zc^^ algorithm [S3] is used as a jet clustering algorithm, thus I 
concentrate on discussing the case of using the anti-fcj. algorithm. The anti-fcj. algorithm basically contains 
two parameters which we can choose their values freely, namely the radius parameter and the lower 
transverse momentum cutoff on jets p’xluf Hence should be defined as 

- y,Y + (</', - </>,)" > , (3.2a) 

Pn > PTc^t) (3.2b) 

where pj’iJ ?/i are the transverse momentum with respect to the beam, rapidity and azimuthal angle 

of outgoing parton i. Imposing cutoffs on the rapidity of partons in not needed. In order to avoid missed 
phase space, the above two parameters must satisfy 

yTlut>PTcut, (3.3) 


(2.14) may 


Let us consider generating an n-jet event according to the n th term in eq. ( |3.1| ), 

cr^x -I-n; {p}x+„ > ^{tx+nAp}x+n ^ {p}X+niets'Ap}x+n) ■ (3.4) 

"^This is the remnant of the fact that there existed the model B too, which turned out not to work properly and thus is not 
described in this paper. 


7 





A parton shower program is executed on an event sample {p}x+n generated by the LO cross section a(X + 
n; {p}x+n > Qraf by setting the shower starting scale to tx+n- Once the shower evolution is performed 
until the shower cutoff scale, all partons in the final state within a rapidity range \y\ < are clustered 
to construct jets according to the anti-fc^ algorithm with the radius parameter and the lower pj, cutoff 


Prcuf If fbe number of the constructed jets is not identical to n, the event sample is vetoed. If the 
event sample survives, the distance parameters Ai? defined in eq. (3.2a) between the n-jet and the n-parton 
are calculated, and then it is checked whether the following is satisfied between each of the n-jet and each 
of the n-parton, or not 


parton < ^ma.tch ^ ■ (3-5) 

If a jet satisfies the above relation with a parton, it means that the jet is close to the parton in terms of 
the measure of the anti-fc^ algorithm. The jet is called matched with the parton. If a matching between a 
jet and a parton is confirmed for all of the n-jet and all of the n-parton, then the event sample is accepted. 
Once the event sample is accepted, all partons in the final state are again clustered to construct jets which 
are defined in the same way as above, but this time only those within a rapidity range \y\ < 
clustered. A value for should reflect actual detectors at experiments and can be different from 

which is used during the merging procedure. The following relation should be satisfied between the two 

(3.6) 

If values for these two cuts are different, the jets constructed at this stage cannot be identical to those 
constructed during the merging procedure and thus are not necessarily matched with the n-parton in some 
events. However, if is large enough so that there is only a small fraction of the events {p}x+n above 

2 /)?®*®®*^, then the effect is expected to be small. In the numerical studies, I use 2/cut'^'^* = 6 and 2/cu’t'* = 7- The 
constructed jets will be used for physics analyses. 

The accepted events have n-jet exclusively and each of the n-jet is close to each of the n-parton produced 
by the LO cross section a{X + n) in terms of the jet measure. Angular correlations between the n-jet should 
follow those between the n-parton. Since angular correlations between the n-parton is at the LO accuracy, 
those between the n-jet should also be at the LO accuracy. 

As is clear, the algorithm calculates the effects of the Sudakov form factors in the same way as the 
MLM [T?lITl] and the fcj^-jet MLM [TTIITB] algorithms]^ The differences from these algorithms are 

• The definition of jets used during a merging procedure is set identical to the one used during physics 
analyses of jets. 

• The LO n,j, 3 ^,j.-parton production cross section is not allowed to produce the events which contain more 
than n,„ 3 ^,,-jet. 

As to the first point, the MLM and the fcj_-jet MLM algorithms use the cone algorithm and the exclusive 
A:j_-jet algorithm, respectively, as a clustering algorithm, and parameters that the clustering algorithm 
contains are chosen independently of the definition of jets used in physics analyses. In other words, the 
definition of jets used in physics analyses is nothing to do with a merging procedure in these merging 
algorithms. In my algorithm, the first point is essential in order to guarantee the LO accuracy of angular 
correlations between jets. 

The maximal number of partons produced by the LO cross section is limited. If we follow the 
treatment of the LO n^^ax'PS'rton production cross section in the MLM and the A:j_-jet MLM algorithms, 
the cross section is allowed to produce the events which contain or more than n,^,^,;-jet and each of the 
hardest njj,j,,;-jet is required to be matched with each of the njj,j,,;-parton. In this approach, the kinematics of 
the additional jets i.e. jets softer than the terms of the jet measure is determined by the leading 

logarithmic (LL) parton shower and hence they do not have angular correlations with any other jets. Since 


®See the discussion at the end of Section!^ 




our objective is to predict angular correlations between jets at the LO accuracy, not only do the jets used 
for an observable such as A;/) = — (j )2 have angular correlations between them at the LO accuracy, but all 

of the jets including jets not used for the observable should have angular correlations between them at the 
same accuracy, because even the jets not used for the observable can affect the observable kinematically. If 
some of the jets not used for the observable are generated by the LL parton shower, then the accuracy of 
the observable will be less than the LO. Considering this argument, the treatment of the LO 
production cross section in the MLM and the fcj_-jet MLM algorithms is not appropriate to our objective. 
In my algorithm, the second point above guarantees the LO accuracy of angular correlations between jets 
in all events. 


The parameter (3.5) is an important parameter, since the Sudakov suppression can depend 

on it. The implementations of the MLM algorithm in Alpgen [ini E] and in MadGraph5_aMC@NLO [34] 


use C. 


match 


= 1.5. However there is nothing that uniquely determines C, 


match ’ 


thus which should be 

considered as a tuning parameter and can also depend on a jet clustering algorithm. In my implementation 
of the new algorithm with the anti-fc^ algorithm, = 2.0 is used. 


Below I describe a procedure of event generation for completeness, 
are performed according to this. 


The numerical studies in Section]^ 


1. Generate the event samples for the X + 0,1,..., rij^jg^'P^rton production processes at proton-proton 


(pp) collisions according to the LO cross sections, i.e. {p}xAp} 


JC + l) ■ 


Ap} 


X+u 


Here the tt 


production is simulated, thus X = tt. MadGraph5_aMG@NLO[35] version 5.2.2.1 is used for this 
purpose. The soft and collinear singularity is regularized by imposing eq. (3.2). I use 

= 20 GeV. (3.7) 


= 0.4, 


A fixed value t^ is used for the scales in the strong couplings and in the parton distribution functions 
(PDFs). The PDF set CTEQ6L1 |36| is used. Note that, when the X + n-jet events are exclusively 
subjects to a study, only the events {p}x+n needed to be generated. 

2. Select an event sample for the X + n partons process, i.e. {p}x+n^ with the probability proportional 
to its integrated LO cross section obtained in step 1, 


a{pp X n) 

X + i)' 


(3.8) 


3. Gonstruct a PYTHIAS parton shower history of the {p}x+n- history consists of intermediate 

events {p}x+(n-i)Ap}x+in- 2 )^''' > {p}x+^^ ’ ’' . Mx+i: {P}x with the clustering scales < 

• • • < t < • • • < ^2 < • 

4. Galculate the weight function for the {p}x+n- define the energy fractions and the parton types 

of the incoming partons in the {p}x+i by and , respectively. The weight function for 

the {p}x+i is given by |20| 


f Ai+iA+iAp} x+i) = 




from which the total weight function for the {p}x+n i® 

n 

Yl f {z,+i,ti+iAp}x+i}, 


(3.9) 


(3.10) 


where = t^^ and The PYTHIAS shower model uses the different strong couplings 

for initial state radiation and final state radiation. Therefore, when a radiation is classified as the 
initial state radiation (the final state radiation) with a scale t^_^_^ by the shower history construction. 
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®s(^z) = 1-37 (1.383) is used for the factor in eq. (3.9). The scale should be determined from the 
intermediate event {p}x- I define it by 


tx = Erp[t) X Erp(t), (3.11) 

where = m? +Py. This scale tx is also used as the renormalization scales of the strong couplings 
for the hard process, namely 


as(^A) 


(3.12) 


is added as a multiplicative factor in the weight function. I use = 0.13 in the above factor. 

Once the weight function is calculated, the {p}x+n i® re-weighted with the function. However, since 
the weight function is not bounded above by unity, the upper bound of the weight function must be 
found at first by calculating the weight function for a large number of events {p}x+n- integrated 
LO cross section obtained in step 1 has to be multiplied by the obtained upper bound of the weight 
function. 


5. Calculate the effects of the Sudakov form factors by using the method described above. I use the parton 
shower model [551 IHl in PYTHIAS [551 HH] version 8186. The default tune of the version 8186, 
tune 4C [26] . is used. The parameters in the anti-fc^^ algorithm and the cutoff values on the rapidity 
of partons are set to 

p^^l, = 30GeY, = y“ = 5. (3.13) 


I use Fastjet m version 3.1.0 for executing the anti-A:.p algorithm. The parton shower starting scale, 
which is the maximal shower evolution scale and is denoted by in 6 < 1 - (3.1 1 , is set to 

tx+n = ^x/4 + ti/16 + ^ 2/16 H-+ ^„/16. (3.14) 


6 . 


Repeat the above procedures from step 2 to step 5 until a large number of the accepted event samples 
are obtained. 


4 Numerical studies 


In this section, the model A algorithm is validated by discussing numerical differences between its predictions 
and the predictions of another merging algorithm. As another merging algorithm, I choose the CKKW-L 
merging algorithm p[ii[5n], which is one of well-established algorithms The top quark pair production 
process at the 14 TeV LHC is simulated, and then angular correlations between the two highest transverse 
momentum jets are studied by using the tt -I- 2-jet events exclusively. The merging scale in the CKKW-L 


algorithm is set as eq. (3.7). 


As an observable measuring angular correlations between the two hardest jets i.e. 
transverse momentum jets, the azimuthal angle difference is chosen 


the two highest 


A(j) = (pi — ^ 


(4.1) 


The events which contain 2-jet are exclusively picked up at first and then if the 2-jet pass the following 
rapidity cut, Ap is evaluated by using the 2 -jet, 


yi X 2/2 < 0, Ivi - 2 / 2 1 > 3.5 or 2.5, (4.2) 

®In my implementation of the CKKW-L algorithm, the method of phase space separation is different from the original one. 
This is the one called the CKKW-L-h algorithm in ref. m- However numerical differences are found small. 


10 







0.010 

0.005 


0.001 
5.x 10"'* 


l.xlO"'* 
5.x 10"^ 




Figure 1: The normalized differential cross section as a function of the of the highest py jet (left panel) and 
that of the second highest jet (right panel) in the 2-jet events. The correspondence between the curves and the 
merging algorithms is labeled inside the left panel. 


which are often called the vector boson fusion (VBF) cuts. One of the 2-jet which has a positive rapidity 
is chosen for (j)i and the other jet which has a negative rapidity is chosen for (^ 2 - The studies based on the 
leading order (LO) ft-h 2-parton production cross section [51 [7] have shown that there are strong correlations 
between the two partons with a large rapidity separation, when the tt is near or far away from the threshold 
of its production |7]. Here I do not impose any cuts on the tt, however. The rapidity cut on jets is set as 
\y\ < 4.5. The tt is assumed stable, since the main purpose of this study is to investigate a way of accurately 
modeling the kinematic activity of jets induced by the hard process. Hence the hadronization after the 
shower evolution and the multiple interaction in PYTHIAS are also turned off. 


At first, I show the transverse momentum pj, distributions of jets. In the left and the right panels of 
Figure the normalized differential cross section as a function of the transverse momentum pj, of the 
highest prp jet and that of the second highest pj, jet in the 2-jet events are shown, respectively. The blue 
solid curve represents the distribution of the model A algorithm. The red solid curve and the black dashed 
curve represent the distributions of the CKKW-L algorithm, the maximal number of partons produced 
by the LO cross section is = 2 for the former and = 3 for the latter. 


The important difference of the two algorithms is that only the 2-parton production cross section 
contributes to the 2-jet events in the model A algorithm, while the 0, 1 and 2-parton production cross 
sections contribute to the 2-jet events in the CKKW-L algorithm at = 2 and the 0, 1, 2 and 3-parton 
production cross sections contribute to the 2-jet events in the CKKW-L algorithm at 


The panels show that the model A algorithm produces softer distributions than the CKKW-L algorithm, 
and that the CKKW-L at = 2 and the CKKW-L at u„jax = 3 produce consistent distributions. Since 
the two algorithms calculate the Sudakov form factors in the very different ways, the differences in the pj^ 
distributions are not unexpected. The model A algorithm contains a tuning parameter in eq. ( |3.5| ). 

In addition to this, the parton shower starting scale, which is defined as in eq. (3.14), is not uniquely 


^The contribution of the 0 and 1-parton cross sections to the the 2-jet events is less than 5%, since the merging scale is set 


as eq. (3.71. 
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Figure 2: The normalized differential cross section as a function of Acj) with the rapidity separation \yi — J/ 2 I > 3.5 
(upper two panels) and |i/i — 1/2I > 2.5 (lower two panels) in the 2-jet events. The correspondence between the curves 
and the merging algorithms are shown inside the panels. 


determined. Therefore, we may tune the algorithm by using these parameters so that the pj^ distributions 
become consistent with the experimental data, before predicting Acj) distributions. This issue is discussed 
again in Section 

In Figure the Acj) distributions are shown. The normalized differential cross section as a function of 
Acj) with the rapidity separation \yi — 2 / 2 1 > 3.5 and that as a function of Acj) with the rapidity separation 
IVi “ 2 / 2 ! > 2.5 are shown in the upper two panels and the lower two panels, respectively. The correspondence 
between the curves and the merging algorithms are shown inside the panels. 

In the left panels of Figure the CKKW-L at = 2 and the CKKW-L at rimax = 3 are compared. 
The difference between the two predictions, particularly in the upper left panel, can be explained from a loss 
of the correlation due to the contribution from the leading logarithmic (LL) parton shower in the CKKW-L 
atn,„ax = 2 M- 

In the right panels of Figure the model A algorithm and the CKKW-L at Umax = 3 are compared. It is 
shown that the model A algorithm produces a consistent prediction with the CKKW-L at = 3, despite 
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that only the LO tt + 2-parton cross section contributes to the events in the model A algorithm. Since the 
LO accuracy of A(j) is guaranteed in the model A algorithm, this observation confirms the finding in ref. [23] 
that the LO accuracy of in the 2-jet events cannot be fully retained due to the contribution from the LL 
parton shower in the CKKW-L at n„iax = 2 and it can be retained by merging the LO ti + 3-parton cross 
section additionally i.e. the CKKW-L at = 3. It should be worth emphasizing that, in the model A 
algorithm, the LO ti -|- 3-parton cross section never contributes the Ac/) distributions in Figure]^ since it 
contributes exclusively to the 3-jet events. 


5 Summary and discussion 

In this paper, a new tree level merging algorithm which guarantees the leading order (LO) accuracy of 
angular correlations between jets is proposed. The new algorithm has been constructed with a goal of 
eliminating the ambiguity in the accuracy of jets, namely it is not necessarily clear whether the kinematics 
of a jet constructed by clustering partons in final state after a merging procedure has the LO accuracy or 
not. This ambiguity potentially exists in merging algorithm, because their virtue is smooth combination of 
LO cross sections and the leading logarithmic (LL) parton shower. The new algorithm requires that n-jet 
events are generated exclusively according to the LO n-parton production cross section and each of the n-jet 
is close to each of the n-parton in terms of the jet measure. As a result, angular correlations between the 
n-jet should follow those between the n-parton and thus the LO accuracy of them is robust. Furthermore, 
as long as the n-jet events are exclusively under consideration, only the LO n-parton production cross 
section is needed and hence event generation is efficient. 

The ti production in proton-proton collisions is simulated and then the distributions of the azimuthal 
angle differences A(j) = — (j )2 between the two highest transverse momentum pj, jets with large rapidity 

separations, namely \yi — 2 / 2 ! > 3.5 and \yi — 2 / 2 ! > 2.5 in addition to x 2/2 < 0 i.e. the VBF cuts, are 
studied by using the 2-jet events exclusively. The new algorithm is evaluated by comparing its predictions 
with the predictions of one of well-established merging algorithms, the CKKW-L algorithm. 

It has been observed that the new algorithm produces a consistent prediction with the CKKW-L at 
rijnax = 3, where denotes the maximal number of partons produced by the LO cross section, despite 
that only the LO ti + 2-parton cross section contributes to the events in the mew algorithm. Since the 
LO accuracy of Acj) is guaranteed in the new algorithm, this observation confirms the previous finding that 
the LO accuracy of A(j) in the 2-jet events cannot be fully retained due to contribution from the leading 
logarithmic (LL) parton shower in the CKKW-L at n^ax = 2 and it can be retained by merging the LO 
ti + 3-parton cross section additionally i.e. in the CKKW-L at = 3. 

In the numerical studies of this work, only the Acj) distributions in the ti + 2-jet events are studied and 
only the CKKW-L algorithm with the PYTHIAS parton shower is used as one of other merging algorithms 
to validate the new algorithm. Therefore, the observation in this work is limited to this situation, obviously. 
However, it is always the case that other merging algorithms including the CKKW-L algorithm potentially 
have the ambiguity in the accuracy of jets, dependently of a process, an observable and a shower model. In 
contrast, the new algorithm does not have the ambiguity and always guarantees the LO accuracy of angular 
correlations between jets, independently of a process, an observable and a shower model. 

There are basically two parameters which cannot be determined uniquely in the new algorithm, namely 
Cmatch definition of the parton shower starting scale. This fact might be seen as a weak point in the 

algorithm. However merging algorithms and parton shower programs are just models after all. A promising 
approach will be to tune the algorithm together with a parton shower model by using these parameters so 
that the pj, and the rapidity distributions of jets in a given process become consistent with the data at first 
and then make the predictions of observables measuring angular correlations between the jets. 
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